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Abstract 

Most methods for modelling dynamics posit just two time scales: a 
fast and a slow scale. But many applications, including many in con- 
tinuum mechanics, possess a wide variety of space-time scales; often 
they possess a continuum of space-time scales. I discuss an approach 
to modelling the discretised dynamics of advection and diffusion with 
rigorous support for changing the resolved spatial grid scale by just 
a factor of two. The mapping of dynamics from a finer grid to a 
coarser grid is then iterated to generate a hierarchy of models across 
a wide range of space-time scales, all with rigorous support across the 
whole hierarchy. This approach empowers us with great flexibility in 
modelling complex dynamics over multiple scales. 



Contents 

1 Introduction 

2 Centre manifold theory supports multiscale models 

3 Coarsen linear advection-dispersion 

3.1 Computer algebra approximates the slow manifold . . 

3.2 The coarse scale, slow dynamics model precisely . . . 

3.3 Coarse scale dynamics are attractive 

3.4 A normal form projects initial conditions 

3.5 Extend elements for a multigrid hierarchy 



10 
11 

12 
13 
17 



* Computational Engineering and Sciences Research Centre, Dept Maths and Comput- 
ing, University of Southern Queensland, Toowoomba, Queensland 4350, Australia. 



1 



4 Multiscale modelling iterates transformations 

4.1 Diffusion emerges from discrete dispersion . . 

4.2 Renormalise advection-dispersion 

4.3 Approximate the nonlinear Burgers' dynamics 



p 

20 
22 
24 



5 Conclusion 



1 Introduction 

Multiscale methods promise efficient computation and simulation of many 
important physical systems [31 e.g.]. Dolbow et al. [10] identify that criti- 
cal applications include fuel cells, subsurface contaminant transport, protein 
folding, climate simulations, and general networks. Most multiscale mathe- 
matical methods for modelling disparate scales presume just two scales: small 
lengths and large lengths; fast times and long times; fast variations and slow 
variations; microscopic and macroscopic pUl EED, e.g.]. Most such methods 
then seek effective models or properties on the large /long/slow macroscales 
by 'averaging/homogenising over' the small/fast microscales [TT| |2"21 e.g.]. 
Here we explore a novel mathematical framework to modelling dynamics 
over many length and time scales; the framework is supported by modern 
dynamical systems theory. 

Of course most two scale modelling methods will also work over many 
scales. The same techniques that construct and support a slow model of 
rapid variations may also apply to construct and support a superslow model 
of the slow variations. The same techniques may then also apply to construct 
and support a megaslow model of the superslow variations; and so on. Cru- 
cially, in most established methods each of these constructions require a large 
'spectral gap'; that is, they require an absolutely clear separation between 
the micro and the macro scales; a parameter such as e measures the scale 
separation, and the requirement for extreme scale separation is provided by 
theorems invoking "as e — > 0". In contrast, multigrid iteration for solving 
linear equations transforms between length scales that are different by (usu- 
ally) a factor of two [51 e.g.]; some variants of multigrid iteration use an 
even smaller ratio of length scales [2E1 e.g.]. Recently Brandt [I] proposed a 
method for molecular dynamics without large scale separation using so-called 
systematic upscaling. Analogously, here we explore modelling dynamics on a 
hierachy of length scales that differ by a factor of two and hence the 'spectral 
gap' is finite and typically much smaller than required by popular extant 
methods for modelling dynamics. Section [2] rigorously supports such models 
with centre manifold theory [HI e.g.]. 
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Recall that multigrid iteration, using restriction and prolongation oper- 
ators, transforms between length scales differing by a factor of two e.g.]. 
The full multigrid iteration involves iterating the restriction and prolonga- 
tion transformations to cross large changes in length scales by taking many 
coarsening or refining transforms where each step changes the resolved length 
scales by a factor of two. Similarly, Section[3]starts our modelling of dynamics 
by exploring a transformation of dynamics from one length scale to another 
that is a factor of two coarser. Section [HI see (fl3l) . demonstrates that to some 
controlled approximation the discrete advection diffusion equation 

= -c±(uj +1 -Uj--,) + d(Uj+i -2uj +uj_i) , (1) 

for evolving grid values Uj(t) on a grid of spacing h (and hence with 'advec- 
tion speed' ch and 'diffusion' dh 2 ) is, on the coarser grid of spacing h. = 2fv, 
justifiably modelled by 



duj 
~dT 

where c = \c and d = |d + , (2) 



-c\ (u j+1 - Uj_| ) + d(u j+1 - 2uj + u M ; 

,2 

16d 



for evolving coarse grid values Uj; these coarse grid values are defined to be 
the fine grid values at every second point on the fine grid, ttj = U2j • In- 
triguingly, the key to the approach is to take one step backwards in order to 
take two steps forward: at any level we embed the dynamics ([I]) in a higher 
dimensional problem, then analysis systematically derives the lower dimen- 
sional, macroscale model The geometric approach to modelling of both 
centre manifold theory [UJ e.g.] and normal form theory [2HJ e.g.] justifies the 
model (EJ) using a similar approach to that of holistic discretisation [26l e.g.]. 
The enhancement of the diffusion by c 2 /(16d) evident in (TSJ) on the coarse 
grid comes from resolving the dynamics on the finer grid in constructing the 
model on the coarser gridQ the enhancement ensures the coarse model ([2]) 
stably models the fine grid dynamics ([T]); intriguingly this enhanced disper- 
sion is precisely that implicit in cyclic reduction, a multigrid method, to find 
an equilibrium of such advection-dispersion problems, but here derived for 
dynamic problems with a different theoretical base. The coarse model (j2J) 
implicitly prescribes a 'restriction operator' that transforms the dynamics of 
advection-dispersion from one grid to another with twice the spacing. 

Others also explore dynamics across space-time scales. Griebel, Oeltz 
& Vassilevski [15] developed space-time multigrid numerics to find optimal 



increasing the dissipation at coarser levels has proven effective in multiscale methods 
for compressible fluid flows p. 10, e.g.]. 
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Figure 1: schematic picture of the multigrid underlying the multiscale de- 
scription of the dynamics: as usual for multigrids, the grids of spacing 
are stacked vertically with relatively coarse grids above fine grids; 
the dynamic variables of the hierarchy of discrete models are uj^. 



control of problems governed by parabolic differential equations. They base 
their cross-scale transformation on an algebraic multigrid. The systematic 
upscaling by Brandt jl] uses multigrid ideas to progressively coarsen atomic 
simulations of polymer folding. These approaches are largely computational 
whereas here we develop algebraic transformations that then are used compu- 
tationally. Another major difference is that the slow manifolds constructed 
here provides a coarsening and interpolation, across length scales, that is 
specifically adapted to the dynamics of the problem rather than being im- 
posed on the problem. Section 13.11 constructs the slow manifolds by system- 
atically approximating exact closures provided by the fine scale dynamics. 

Section H] explores iterating our transformation to model dynamics across 
each and every intervening length scale. For example, repeating the transfor- 
mation from fine ([T]) to coarse ([2]) gives a hierarchy of models all of the form 
of the advection-dispersion equation (CQ) but with differing coefficients. At 
the £th level, with grid spacing h' £ ' = 2 { h, the corresponding grid values ui 
evolve according to (ED) but with coefficients c^' and d' £ ' determined by the 
recurrence 

c^ = ±c^ and ^ +1 » = X' + ^. (3) 
On successively coarser grids the coefficients thus are 

c w = 4 and d ( « = J£U< fl where d (£+1) = \ ( d w + J-V (4) 
2« 2 f+1 2 V dW / 
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Observe that d' £ ' — > 1 quickly as I increases (as (j3J) is equivalent to Newton's 
iteration to find the zeros of d 2 — 1); hence as the grid coarsens, the £th level 
model quickly becomes simply the upwind model 



du 




■j-i 



dt 



— c 



when c > . 



Our multigrid modelling transformation naturally recognises that advection 
dominates diffusion on coarse grids: the cross scale transformation, the map 
from fine ([T]) to coarse ([2D as summarised by , not only preserves the advec- 
tion speed, but also models the advection in a stable scheme that preserves 
non-negativity. Further, in the absence of advection, c = , the transfor- 
mation (J3j) preserves the effective diffusion across all scales: d' £+1 ' = . 
These are some simple results. Section [H explores further issues in transform- 
ing both linear and nonlinear discrete dynamics across many scales. 

The centre manifold and normal form [121 El 1201 e.g.] approach estab- 
lished here provides a framework for dynamical modelling that links what are 
conventionally called multigrid [SI e.g.], wavelets 0, e.g.], multiple scales [231 
e.g.], and singular perturbations [361 e.g.]. This framework applies to not 
only the linear dynamical systems that are the main focus of this article, 
but also applies to nonlinear systems [2S1 CHI e.g.] and to stochastic sys- 
tems [3 HI EO, EH e.g.]. Here, because it is simplest, we focus on transforming 
dynamics within the same algebraic form, but in principle the methodology 
can support the emergence, via nonlinear interactions, of qualitatively differ- 
ent dynamics on macroscales (as promoted by the heterogeneous multiscale 
method [TTl e.g.]). By rationally transforming across both space and time 
scales, a long term aim of this approach is to empower efficient simulation 
and analysis of multiscale systems at whatever level of detail is required and 
to a controllable error. 

This approach to transformation from one scale to anther may in the 
future illuminate complex systems simulations on both lattices and with cel- 
lular automata. 

2 Centre manifold theory supports multiscale 
models 

This section establishes new theoretical support for coarsening dynamics from 
a fine grid to a coarse grid of twice the spacing. Suppose the fine grid has 
grid points Xj, spacing h as shown in Figure [21 and has grid values Uj(t) 
evolving in time. The figure also shows the coarse grid points Xj = x 2 j , 
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Figure 2: schematic picture of the equi-spaced fine grid, Xj, with spacing H, 
and grid values Uj . The coarse grid, Xj, with spacing fx = 2h, and grid 
values Uj is superposed. 



spacing h = 2h , and the definition of the evolving coarse grid valued 

u,(t) =u 2j (t). (5) 

Mostly, an overbar denotes variables and operators on the coarser grid, and 
unadorned variables are those on the finer grid. Using overdots to denote time 
derivatives, the question is: how do we transform the evolution ii-j = £uj , 
for some fine scale local operator C, to a coarse evolution Uj = Cuj on the 
coarse grid? 

The theoretical support for multiscale modelling outlined by this section 
applies equally well to nonlinear dynamics: Section 14.31 briefly explores the 
specific nonlinear advection-dispersion of a discrete Burgers' equation. As- 
sume the fine spatial grid is periodic with ru grid points: that is, for definite- 
ness assume the grid is periodic in space x with period rah. For conciseness, 
write equations in terms of centred mean (x and difference 6 operators [2"Tl 
Ch. 7, e.g.] acting on the fine grid. Thus the advection-dispersion equa- 
tion (CQ), but now including some 'nonlinearity' fj with some parameters e, 
is 

Uj = { — c\xb + dS 2 }uj + fj(u, e) . (6) 

I give three illustrative examples of such nonlinearity: a local reaction could 
be prescribed by f j = euj — u? ; a nonlinear advection by f j = Uj|x5uj/(2h\) as 
in the discrete Burgers' equation ( 14T|) : whereas linear diffusion in a random 
medium could be encompassed by fj = e6(Kj6uj)/H 2 for some stochastic 
diffusivities Kj. 

2 This is the amplitude condition that ensures the relation between the coarse grid values 
and the fine grid values are unique and well defined. Others, such as Brandt [3], prefer 
using averages of the fine grid values as coarse variables which would also be acceptable. 
Centre manifold theory allows the coarse variables to be any reasonable measure of the 
amplitude of the fine grid variables. 
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Figure 3: schematic picture of the equi-spaced fine grid, Xj, with spacing H, 
and the coarse grid, Xj, with spacing h = 2h. Three finite elements of 
the coarse grid are shown in exploded view to illustrate their overlap with 
neighbouring elements. Within each finite element new variables Vj t(t), —2 < 
i < 2, replace the fine grid variables u,-(t). 



Centre manifold support We now describe how to support and construct 
the model on the coarse grid of the fine scale, nonlinearly modified, advection- 
dispersion dynamics. 

Analogous to holistic discretisation of pdes [2ol [T5| e.g.], divide the m- 
periodic fine grid into m/2 overlapping elements. Notionally let the jth coarse 
element stretch from Xzj-z to xzj+2 as shown 'exploded' in Figure El As 
shown, denote the evolving fine grid field in the jth element as the 5-tuple 
Vj = (vj _2, Vj _i, v^o, v^i, Vj t 2) , so that at this stage we have just renamed 
the fine grid variables, Uzj+i = v^. Note that the elements overlap: the 
fine grid values U2j±i = Vj±i = Vj±i )=F i ; this overlap empowers us to cou- 
ple the dynamics in neighbouring elements to derive consistent models as 
similarly derived for holistic discretisation [29]Jf| The interelement coupling 
conditions (jHJ) determine the fine grid values Vj ±2, at the extremes of each 
element, and so these are not extra dynamic variables. But, importantly, 
consider the overlapping fine grid values Vj ±i and Vj±i jT i as independent 
dynamic variables satisfying the fine scale discrete equation ©, namely 

v jfi = {-c^6 + d5 2 }v j)i + f j (v j) e), i = 0,±1, (7) 

where these differences and means operate over the fine grid index i. In 
essence I extend the dynamics of the m fine grid variables Uj(t) by an extra 

3 The overlapping elements may be analogous to the 'border regions' of the hetero- 
geneous multiscalc method [TTJ e.g.] and to the 'buffers' of the gap-tooth scheme [551 
e-g.]- 
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m/2 variables. This is the 'one step backwards' referred to in the Introduc- 
tion: in order to rigorously support the modelling of the m-dimensional fine 
scale dynamics by m/2 coarse scale variables, I embed the fine scale system 
in the 3m/2-dimensional dynamics of these overlapping elements. Section [33] 
shows how to choose these m/2 extra degrees of freedom to make forecasts 
from any given fine grid scale initial condition. 

Also analogous to holistic discretisation of pdes [2EJ HEJ e.g.], couple 
neighbouring elements with the conditions 

Vj,±2 =?Vj±i,o+ (1 -?Ho, (8) 

where the coupling parameter y controls the interaction and information flow 
between elements: 

• when y = 1 the elements are fully coupled and the condition ([8]) reduces 
to the statement that the extrapolation of the jth element field to the 
neighbouring coarse grid points, Vj ±2, is identical to the neighbouring 
coarse grid values Vj±i t o{= Uj-ti); 

• when y = the elements are completely isolated from each other and 
thus, linearly, the new fine grid values Vj t evolve quickly to be constant 
in each element. 

This equilibrium when y = , or space of equilibria depending upon the non- 
linearity fj, forms the base for the slow manifold model which when evaluated 
at y — 1 gives the desired model for the fully coupled dynamics. 

I use the overbar in y because it moderates information flow between 
the elements forming the coarse grid. By working to an error (9(y n ) we 
account for interactions between the dynamics in an element and its n — 1 
neighbours on either side. Thus we transform local dynamics on a fine grid 
to local dynamics on a coarse grid as in other multiscale approaches jU e.g.]. 
The size of the locality depends upon the order of error in the coupling 
parameter y. 



The decoupled dynamics have a useful spectral gap Set y — to 

decouple the elements, and neglect the nonlinearity by linearisation. Then, 
independently of all other elements, the linear dynamics in the jth element 
are governed by the differential-algebraic system 
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Seeking solutions proportional to e At this set of linear odes has three eigen- 
values and three corresponding eigenvectors: 



A = 0,-2d,-4d; and 
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(10) 



From these, any zig-zag structures within an element decay exponentially 
quickly, and hence these decoupled dynamics results in constant solutions in 
each element arising on a time scale of 1/d. Over all the m/2 decoupled 
elements these piecewise constant solutions form an m/2 dimensional linear 
subspace of equilibria, the so-called slow subspace, in the 3m/2 dimensional 
state space of the fine grid values Vj and Vj ±1 . Centre manifold theory for 
deterministic systems [HI EQl e.g.] or for stochastic systems (2j (TJ §8.4, e.g.] 
then assures us of the following three part theorem. 

Theorem 1 For some domain of finite non-zero coupling parameter y , and 
if nonlinear, some neighbourhood of the origin in (u, e): 

1. there exists an m/2- dimensional, invariant slow manifold M. of the 
coupled dynamics of the discrete nonlinearly perturbed, advection- dis- 
persion ((7j) with coupling conditions ([8]) — with a dimension correspond- 
ing to each of the m/2 coarse grid elements; 

2. the dynamics on the slow manifold Ai are approached exponentially 
quickly, roughly like exp(— 2dt), by all initial conditions Vj ; i(0) of the 
fine grid values in some finite neighbourhood of M. — that is, the slow 
manifold dynamics faithfully model for long times generic solutions of 
the coupled dynamics; 

3. we may construct the slow manifold model to some order of error 
in y, |e| and |u| by solving the governing, nolinear, discrete advection- 
dispersion (Ej) with coupling conditions (EJ) to residuals of the same 
order. 

Two broad cases arise: if the nonlinearity f j = whenever Vj ; i is in- 
dependently constant in each element — for example the Burgers '-like 
nonlinearity fj = Uj|x5uj/(2h.) — then the approximation is global in 
the coarse grid variables u; alternatively, whenever f j ^ for Vj t in- 
dependently constant in each element — for example the reaction fj — 
euj — uf — then the approximation is local to the origin in u 
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Like systematic upscaling [H pp. 6, 9] and other multiscale methods, this 
approach uses equilibrium concepts. But one crucial difference is that centre 
manifold theory guarantees that the same separation of dynamics occurs in 
a finite neighbourhood about equilibria and hence supports the separation 
of coarse scale dynamics from the fine scale occurs for nontrivial dynamics. 
This approach provides a systematic alternative to the heuristic Fourier or 
wavelet decompositions for a 'local mode analysis' [H §8]: here the local 
modes are determined by the the dynamical system itself through the shape 
of the slow manifold. 

Finite domain After constructing an approximate slow manifold model, 
we evaluate it for coupling parameter y = 1 to recover a coarse grid model for 
the fully coupled dynamics on the fine grid. Is y — 1 in the 'finite neighbour- 
hood' of theoretical support? It is for the analogous holistic discretisation of 
the Burgers' pde [2E]. Similarly, Section |3~B1 demonstrates that the fully cou- 
pled case, y — 1 , is indeed within the neighbourhood of theoretical support 
for the linear (fj = 0) dynamics of flZ|. 

3 Coarsen linear advection-dispersion 

Using the theoretical support of centre manifold theory established by the 
previous section, this section analyses linear advection-dispersion to provide 
the multiscale modelling results summarised in the Introduction. 

3.1 Computer algebra approximates the slow manifold 

Elementary algebra readily constructs general slow manifold models [251 [TBI 
e.g.]. We solve the fine grid, linear, discrete, advection-dispersion equa- 
tion ([7j) with coupling conditions (jBJ) by seeking solutions parametrised by 
the evolving coarse grid values Uj(t): 

Vj >t (t) = Vj^u.y) such that ttj = Gj(u,-y) , (11) 

for some functions Vj ( i and Gj to be determined by the iterative algorithm [32J. 
The base approximation is the slow subspace of equilibria: 

Vj ) i(t) = Vj^u, y) ps Uj such that Uj = Gj(iij,y) ps . 

Computer algebra code [32] systematically refine these slow manifold approx- 
imations. The refining iteration is based upon the residuals of the discrete 
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equation (|7|) with coupling conditions (jSJ). Centre manifold theory then as- 
sures us that the error in approximating the slow manifold model is of the 
same order in coupling parameter y as any remaining residual. For example, 
to errors 0(y 2 ), computer algebra [32] constructs the slow manifold A4 in 
the jth element as 

y) +7E- 1 



(1 
1 - 
1 

1 + +_ 
.(1 -y)+?E 



y(l 



c 

4d? 
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4d? 



Ui 



o(y 2 ), 



(12) 



in terms of the coarse grid centred difference and mean operators, 6 and p.6, 
and shift operator E (define Euj = Uj + i or equivalently E = E ). The terms 
in (|12[) which are independent of advection c, for the fully coupled y = 1 , are 
classic quadratic interpolation from the surrounding coarse grid values u,. 
The terms involving advection, flagged by c, arise through accounting for 
the dynamics of the fine grid values U2j±i and their interaction with the 
surrounding grid values. Equation ( fl~2l) corresponds to the multigrid pro- 
longation operator, but here it is derived by accounting for the fine scale 
dynamics rather than being imposed. 

The evolution on the slow manifold (TT2l) is then the coarse grid model 



y 



-\c\ib + 



^ d + T6d 



u j + o(y 2 ) 



(13) 



Neglecting the C(y 2 ) error, evaluate (fl3"|) at the physically relevant coupling 
y = 1 to deduce the coarse grid model ([2D discussed in the Introduction. 



3.2 The coarse scale, slow dynamics model precisely 



Consider further the linear advection-dispersion (I7j) with interelement cou- 
pling conditions ([8]). The previous subsection constructed an approximation 
to errors O (y 2 ) ; this subsection gives exact formula for all coupling y. 

Seek solutions with structure within the finite elements of the formal op- 
erator form Vj = exp(t/C TL )e n , where, generalising ( |T0l) to non-zero coupling, 
/C n is the nth 'operator eigenvalue' of the advection-dispersion ©-([H]) and 
e n is the corresponding 'operator eigenvector'. Elementary algebra for any 
coupling y reveals the three operator eigenvalues are precisely 



/C li3 = 2d^-1± Wl+y 



1 + 



4d 2 



2d 



u5 



and JCz 



-2d. 

(14) 
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The smallest (least negative) of these, namely /Ci, governs the longest time 
scales in the coupled dynamics. For example, the Taylor expansion^ in the 
coupling y of operator /Ci (the plus case above), upon using the identity 
p 2 = 1 + \b 2 , agrees with the C(y 2 ) evolution ( fl3l) . to the C(y 3 ) approx- 
imation ( l40i) . and to the 0(y 5 ) approximation ( 13T1) of isotropic dynamics. 
That is, the coarse grid evolution operator £ = JC-\ . 

Lemma 2 The coarse grid operator £ = /Ci reproduces exactly the fine grid 
operator of linear advection- dispersion when the elements are fully coupled. 

Proof: In the fully coupled limit, y = 1 , the three operator eigenval- 
ues ( !T4l) reduce to 

L=K\ = d(2fi-2)-|c6, /C 2 = -2d, /C 3 = -d(2p. + 2) + ±c6 . (15) 

Relate to the fine grid operators, via the coarse grid shift operator E and the 
fine grid shift operator E( = E 1 / 2 ), by observing 

2fi - 2 = E 1/2 + TT 1/2 — 2 = E + E -1 — 2 = S 2 
and ±b = \(V/ 1 -t-V 1 ) =1(E-E- 1 ) = ji6 . 

Consequently, the coarse grid model 

tlj = £tij = [d(2p. — 2) — -^cSjuj = [d6 2 — C(j.6]uj , 

is precisely the fine grid, linear, advection-dispersion equation ©, except 
that, having half the grid points, it does not resolve the fine scale, high 
wavenumber, spatial structures that the fine grid can resolve. 6 

Thus the operator £ does indeed model on the coarse grid all the coarse 
dynamics inherent in the fine grid advection-dispersion dynamics. It is only 
the approximation of £ by a truncated Taylor series, such as in the 0(y 2 ) 
model ( TT3l) . that induces errors in the coarse scale model of the long term 
dynamics of linear advection-dispersion. 

3.3 Coarse scale dynamics are attractive 

Consider the spectrum of the advection-dispersion dynamics implicitly de- 
scribed by the operator eigenvalues (fT4|) . On any regularly spaced grid, the 
centred mean and difference operators act on Fourier modes as 

^ ikj = cos(k/2)e ikj and 6e ikj = 2isin(k/2)e ikj 

4 The Taylor expansion converges provided the bracketed coefficient of the coupling y 
is of magnitude less than one. This convergence occurs provided the differences p.6 and 6 
are small enough; that is, for solutions varying slowly enough across the grids. 
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for a component of spatial wavenumber k relative to the grid spacing; the 
wavenumber domain is —n < k < n. Upon taking the discrete Fourier 
transform, the operators p. and 5 thus transform to cos(k/2) = cosk and 
2isin(k/2) = 2isink, respectively, as the fine grid wavenumber k = k/2 
in terms of the coarse grid wavenumber k. Thus from (|T^|) the advection- 
dispersion dynamics on the fine grid elements has spectrum 



M,3 



2d -1± Jl+y 



1 + 



4d 2 



7 V • u 'i • 1<- V 

sm 4 — T- — 2 sm 4 cos 4 
2 2d 2 2 



and A2 = —2d , 



(16) 



for coarse grid wavenumbers |k| < n. Extensive numerical computations 
strongly suggest that < D^^A < 1 , where \f- denotes the square root 
in ( 1161) . for all wavenumbers |k| < 7t, for all coupling < y < 1 , and for 
all advection relative to diffusion, c/d. Consequently, the numerics suggest 
the spectral ordering < SPCK^ < ^CK\ < is maintained across the 

whole relevant parameter domain. Thus, not only does the coarse grid model 
\ij — C\Lj = /C-|Uj accurately model the fine grid dynamics, the coarse grid 
model is the slowest dynamics of the fine grid advection-dispersion. 

Theorem^ ensures an atttractive slow manifold exists in some neighbour- 
hood of coupling y = ; the spectrum ffTUl) demonstrates that the neighbour- 
hood extends to include the case of fully coupled elements, y = 1 . 

We usually cannot construct slow manifolds exactly, as done above; in- 
stead we usually approximate slow manifold by a multivariate power series. 
Thus the practical issue is not just whether a slow manifold exists, but how 
well a truncated power series approximates the slow manifold. Elementary 
algebra shows that a Taylor series of fll6p in y converges at y = 1 provided 



[0 



sin 



: k + 4£ 2 ] sin 2 k < 1 , 



where <£ — c/(2d) measures the advection relative to the dispersion. For all 
parameter £ there is a finite range of small wavenumbers k satisfying this 
inequality. This argument leads to the following lemma. 

Lemma 3 Finite truncations of the Taylor series of the slow operator eigen- 
value /C1 provide accurate approximations of the evolution of the coarse grid 
variables provided the solutions vary slowly enough across the grid. 



3.4 A normal form projects initial conditions 

Suppose we know the fine grid values Uj(0) at the initial time t = 0. This 
subsection addresses the question: what coarse grid values should we give 
to Uj(0) for the coarse grid model to make accurate long term predictions? 
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The obvious answer is wrong [2H El e.g.]: even though we define 
Uj(t) = U2j(t) , we nonetheless should not set the initial Uj(0) = U2j(0). 
The reason is that the transient dynamics of the subgrid scale dynamics 
modifies the appropriate initial value for Uj(0); this modification is sometimes 
called 'initial slip' in physics [T4" l [T3 l e.g.]. In this subsection, a normal form 
coordinate transform of the fine grid dynamics clearly displays the correct 
initial conditions for the coarse dynamics. 

In this discussion, restrict attention on initial conditions to the fully cou- 
pled case of coupling y = 1 . This restriction simplifies by avoiding the 
complicating detail of having variable y, and it focusses on the physically 
relevant case of full interelement coupling. 

Consider the spectral decomposition of the dynamics of the fine grid of 
all the elements. Transform the fine grid evolution to its coarse grid 'normal 
form' of the spectral decomposition 



v,(t) 



ej,iUj(t) 



where vu, = /C]U 



- ej, 2 Vj(t) + e 3 - )3 Wi(t) , 



for the operators /C n in (j!5p and for intraelement structure operators 



\c- d" 


\c + d 



and 



c(E- | -1) + d(E- 1 + 1 

/C n + 2d 
±c(E-1) + d(E + 1) 



(17) 



for n = 1 , 3 . I do not record the two extreme components ±2 in these ej >n 
as Vj t ±2 are identical to Vj±i ( o when fully coupled, y = 1 . Within each of the 
fully coupled elements, a formal expression for the complete evolution on the 
fine grid is thus 



v, t 



e 5>1 exp(t/Ci)uj(0) + e j)2 exp(t/C 2 )vj(0) + e j)3 exp(t/C 3 )Wj(0) , (18) 



for some constants Uj(0), Vj(0) and Wj(0). For example, from (1151) . when 
advection c = the intraelement structure operators simplify to 













" 1 " 


e j)2 oc 





and e^ n oc 


±p- 




±1 




1 








1 



where this last approximate equality holds for fields varying slowly enough 
along the grids. Thus e^j ~ (1,1,1) represents the smoothest variations 
within each element, whereas e^z ~ ( — 1 , 0, 1 ) and e 3 - )3 ~ (1,-1,1) represents 
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fine grid scale fluctuations within an element Jf] Since these fine grid scale 
fluctuations decay rapidly in time t, the long term slow dynamics on the 
slow manifold is just the restriction of ( TT8]) to v, = Wj = , namely 



Vj(t) 



e^j exp(t/Cit)uj(0) 



(19) 



We must choose the initial condition, Uj(0), for the coarse grid values so that 
this evolution exponentially quickly equals the fine grid dynamics Uj(t) = 
exp(t£)uj(0) from the specified initial condition. Such a choice for the initial 
coarse grid value Uj(0) then realises the theoretical promise by Theorem [1]- 
[2] of long term fidelity between coarse grid model and fine grid dynamics. 
Elementary linear algebra determines the coarse grid values Uj(0) through 
evaluating the general solution (ITS]) at time t = , 



Vj(0) = e^iljtO) + e j)2 Vj(0) + e j3 w,(0) , 
and then take the inner product with the left eigenvector 



(20) 



to deduce the following lemma. 



jc+d " 
/Ci +2d 
-\c + d 



Lemma 4 For linear advection- dispersion, the initial coarse grid values are 

zj.i • Vj(0) 



in terms of specified fine element values v,(0). 



(21) 



Despite the definition that the coarse grid values Uj(t) = U2j(t) , the normal 
form coordinate transform accounts for dynamics in fast time initial tran- 
sients so that the correct initial conditions for the coarse grid model is the 
nonlocal and weighted projection (12"TT) . 

The initial condition mapping ( 12~TI) relates to multigrid iteration. When 
advection c = 





" 1 " 




"1" 


Zj,i oc 


2ji 




2 




1 




1 



5 In essence, the element eigenvectorsej in are the natural wavelets [SI e.g.] for the 
advection dispersion dynamics. The difference here is that, being adapted to the dynam- 
ics, the precise shape of the element eigenvectors depends upon the coupling with the 
neighbouring elements. 
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and in the case of slowly varying grid values, this projection from the fine 
grid initial values Vj(0) to the coarse grid initial values is the classic multigrid 
restriction operator [5j e.g.]: namely, that the coarse grid value is the average 
of the nearest fine grid values with a weighting of 1 : 2 : 1 . 



Uniquely prescribe fine element values We have an additional com- 
plication: on the fine grid, the odd grid values U2j±i are shared between 
two neighbouring elements. The grid value U2j+i is represented as both 
Vj i and Vj+i _i, and both of these variables are treated as separate inde- 
pendent variables in the dynamics on each element. We must resolve this 
separation. 

Two independent suggestions resolve the separation with the same result. 
My first suggestion to avoid conflict between the values of Vjj and Vj + i _| is 
to require that Vj i = Vj+i _| at the initial time. The shift operators rewrite 
this identity as 

Ev J|0 = EE- 1 v j ,o, (22) 

where the coarse grid shift E operates on the coarse grid, first subscript 
of Vj^, whereas the fine grid shift operator E operates on the fine grid, 
second subscript of Vj t. For a domain with m fine grid points, that is, 
m/2 coarse grid elements, the compatibility condition (1221) provides an addi- 
tional m/2 constraints to determine uniquely the 3m/2 initial values Vj^fO) 
within the elements from the m fine grid values u-j(O). My second sugges- 
tion is to choose Vj^i (0) so that the unphysical intermediate mode vanishes 
in the solution (fTBl) . that is, so that Vj(0) = in the solution (ITS]) . Then 
there will be no intermediate scale dynamics exp(— 2dt) and the approach 
to the slow manifold model will be the quickest: the only rapidly decaying 
mode will be the ej i 3exp(/C3"t)'vVj(0) mode which, from the spectrum ( 1161) . 
decays more rapidly than exp(— 2dt). Now relate Vj(0) directly to Vj(0) by 
multiplying (12U1) by the left eigenvector corresponding to K,z namely 



Zi,2 = 



L2 



±c(E-1)-d(£ + 1) 


Jc^-IJ + dfE^+l) 



Thus, noting Vj^i = E ±1 Vj o , 

Vj(0) ex {[ic(E-1)-d(E + 1)]E- 1 + [Ic(E- 1 -1) + d(E- 1 + l)]E}v ji0 
= { [\c - d) (EE- 1 - E) + ( \c + d) (E- 1 E - E- 1 ) } v,, . 

Consequently, ensure the mode exp(— 2dt) does not appear at all, Vj(0) = , 
by requiring (EE -1 — E)Vj,o = which is precisely ( 1221) . and by requiring 
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(E E — E )vjo = which is again f[2"2l but just shifted to the left by the 
multiplication by the coarse grid shift E -1 . Thus the condition (1221) ensures 
that neighbouring elements agree at their common points and that the slow 
manifold, long term model is approached quickest. 

Lemma 5 Choosing the embedding to 3m/2- dimensions to satisfy ( 1221) at 
the initial time ensures that (122|) is satisfied for all time in the linear, advection- 
dispersion dynamics on the fully coupled finite elements. 

3.5 Extend elements for a multigrid hierarchy 

As discussed briefly in the Introduction, we aim to transform dynamics across 
a wide range of space-time scales using the multigrid hierarchy illustrated in 
Figure HJ The Introduction used a model of O (y 2 ) , see §3.11 to transform 
advection-dispersion on a fine grid to advection-dispersion of the same form 
on a coarser grid. This transform iterates simply across all scales. How- 
ever, when we seek more accuracy, say errors O (y n ) for n > 2 , the lin- 
ear advection-dispersion dynamics transforms into a model of the form 
dtlj/ dt = G(Uj_ n+ i , . . . , Uj +n _i , y) that involves 2(n— 1 ) neighbouring coarse 
grid values. For example, to errors 0(y 4 ), fine scale isotropic dispersion 
(equation ([6]) with c = f j = 0) transforms to the coarser scale dispersion [32] 



that through 6 4 and 6 6 involves Uj±2 and Uj±3 . Consequently, to empower 
us to transform coarse models over a hierarchy of grids we must widen the 
elements defined in Figure E] to include more fine grid points. This subsection 
widens the elements while maintaining the spectrum (I1UI) ensuring the centre 
manifold support P, [20j e.g.]. 

This subsection, as seen in equation (1231) . avoids the overdots for time 
derivatives as we invoke different time scales on each level of the hierarchy. 
Interestingly, it eventuates that not only do we overlap the elements, but 
also, in some sense, overlap the time scales. 

The general form of linear dynamics on a grid Suppose at some 
level of the multigrid hierarchy we know the discrete operator governing the 
evolution of grid values u,-(t). Decompose the discrete evolution operator as 
the sum 



duj 



(23) 



dt 



duj 



[£, + C 2 + £ 3 + • • • + Ai-l] Uh 



(24) 



dt 
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where the kth discrete operator Cy_ has stencil width 2k+1 ; that is, C^Uj only 
involves Uj-k, . . . ,Uj + ^. This decomposition terminates, as written in 024|) . 
when we restrict attention, by working to errors 0(y n ~) , to operators of some 
maximum finite width. The decomposition is not unique as specified so far; 
however, as apparent in ( |23i) . a specific unique decomposition naturally arises 
when we generate the models to errors 0(y TL ^ in some coupling parameter y. 
Thus suppose there is a natural 'ordering' parameter y such that, instead 
of (124j) . the discrete evolution equation may be written 

^ = [A + yC 2 + y 2 C 3 + ■■■+ y n - 2 £n-!] u } . (25) 

At all levels, except the very finest level 0, this natural parameter y is to 
be the coupling parameter of the finite elements of the grid one level finer 
than than current level. As always, we suppose that evaluation of ( 125]) at 
y = 1 gives the physically relevant model ( l24l) . whereas y = provides 
a base for theory to support models at non-zero y. Additionally insisting 
on the operator L\ being conservative implies £i must represent advection- 
dispersion dynamics and implies that the y = dynamics, duj/dt = £iUj , 
provides the same sound base for applying centre manifold theory. Note that 
the coupling parameter of finite elements at the current level is still y. That 
is, still couple neighbouring elements with the condition (jSJ). 

Anticipating the support by centre manifold theory, derived in a couple of 
paragraphs, we expect to construct a coarse grid model of the dynamics 025p 
in the form 

^ = [yC, + y 2 C 2 + y 3 £ 3 + ■ ■ • +y n ~ 1 £n- 1 ] Uj + 0{y n ,y n -') , (26) 

for some coarse grid operators (implicitly a function of the artificial y) 
which will be of stencil width 2k + 1 as the parameter y counts the number 
of interelement communications. The renormalising transformation requires 
two extra ingredients: first remove the fine grid ordering by setting y = 1 (so 
operators Cy_ are no longer a function of y); and second introduce a coarse 
grid time scale t = t/y (which is the same time when y = 1), then, upon 
dividing by y, the coarse grid dynamics become 

^ = [A + yC 2 + y 2 C 3 + ■■■+ y n - 2 £n-i] Uj + 0{y n - ] ) . (27) 

The coarse model (127|) has exactly the same form as the fine model ( 1251) . 
By introducing the coupling (jSJ) across all levels of the hierarchy, and by 
introducing a hierarchy of times, which all collapse to the same real time 
when y = 1 , and working to some order of error in coupling, models of 
the form (125j) are transformed and renormalised across the entire multigrid 
hierarchy. 
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Widen the elements Assume we wish to construct slow manifolds to 
errors 0(y n ,y n_1 ) with the aim of using centre manifold theory to support 
the modelling of ( |25|) by ( l2?j) . Extend Section [2] by widening the jth element 
to extend over the interval [xj_ n , Xj +n ] and also to possess the (2n + 1 ) fine 
grid variables Vj = (vj _ n , . . . , Vj >n ) . These extra variables are not extra 
degrees of freedom. Let these fine grid variables evolve according to 

dv^ ^ -> ^ 

= /:iv j , i + y£ 2 v j , i + y 2 £ 3 v j>i + ---+y n - 2 /: n _ 1 v j , i , |i|<n, (28) 

where we adopt the unusual convention that when applied within the ele- 
ments, the operator L\y^\ is its original definition when |i| + k < n but 
is zero otherwise (for |i| + k > n). Adopting this convention ensures that 
the operators on the right-hand side of (|28p do not 'poke outside' of the 
jth element; in effect, this convention truncates the sum in (|28[) to remain 
within the jth element. Such truncation incurs an error 
evolution of a variable v^. However, as variable v^i only affects the crucial 
central core variables of the element, v^v for < 2 as shown in Figure [3], 
via terms of (^(y' 1 ' -2 ), the net effect of this conventional truncation is an er- 
ror C?(y n_1 ) which is the same as the assumed order of error of the analysis. 
The coupling condition (jHJ) closes the dynamics on these widened elements. 
In essence we do not have new dynamics outside of the central core of each 
element, instead, in effect, we simply extrapolate the dynamics to the outside 
of the central core. 



Centre manifold theory support When the fine grid 'ordering' param- 
eter y = and interelement coupling parameter y = the dynamics on the 
m/2 elements reduces to 

— £i v j,i f° r |t| < ft , and Vj j± 2 = Vj >0 . (29) 

Each element is decoupled from the others. The general conservative, linear, 
operator is the advection-dispersion operator, C-\ = — cp.6 + dS 2 for some 
constants c and d. As for the earlier ([9]), for each of the extended elements 
there are still precisely three eigenvalues of (}2"9"]) . namely A = 0,— 2d,— 4d. 
Corresponding eigenvectors are the constant Vjt oc 1 , the artificial v, ^ oc 
sin(i7r/2) [(1 +<£)/( 1 — <T)] t/2 , and the zig-zag mode Vj,i oc (— l) 1 . Conse- 
quently, centre manifold theory implies Theorem [1] also applies to the sys- 
tem (1281) with coupling conditions (JHJ) to ensure: firstly, that an m/2 di- 
mensional slow manifold exists for the dynamics of the coupled elements; 
secondly, the coarse scale dynamics on the slow manifold are attractive; and 
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thirdly, that we may construct the slow manifold to any desired error — this 
section assumes errors O (y n ~^ , y n ) . 

The next subsection proceeds to briefly explore the resultant models of 
advection and dispersion over a hierarchy of multiscale grids as supported by 
this theory. 



4 Multiscale modelling iterates transforma- 
tions 

This section explores three example applications of transforming dynamics 
repeatedly across the wide range of length and time scales on a multigrid hi- 
erarchy. Section |4~T1 shows how continuum diffusion emerges from microscale 
dispersion. Section 14.21 deonstrates that the nonlinear transformation from 
one scale to another of linear advection-dispersion has a fixed point of a sta- 
ble upwind model. Section 14.31 discusses briefly the transformation of the 
nonlinear Burgers' pde. 



4.1 Diffusion emerges from discrete dispersion 

The multiscale modelling of discrete dispersion, when the advection coeffi- 
cient c = , reduces to a remarkably simple linear transformation. Here 
we explore the exact slow manifold transformation from a fine grid to a 
coarser grid. Iterating this transformation proves that, in the absence of ad- 
vection, the continuum diffusion equation naturally emerges very quickly on 
macroscales. 

Linear dynamics which are left-right symmetric (isotropic) can be ex- 
pressed in terms of only even order central differences. Our slow manifold, 
multiscale modelling preserves this form. Suppose the evolution at grid level I 
is governed by 



duf } 



^y^cMs^uf 1 , (30) 



=1 



for some coefficients c p ; for example, the second order difference coefficient 
Ci = d used earlier. In practical constructions, invoking an error d}(y n_1 ) 
truncates to a finite sum this 'in principle' infinite sum. Computer alge- 



bra [32] , supported by the theory of Section 13.51 derives the dynamics at the 
next coarser level of the multiscale hierarchy, namely 



du) 



20 



y5^2 t ^ i28 c 2 + 64 c 3 y° 



+ y 3 



5 -W , 5 JD 
1 6384^1 ' 4096 ^2 



3 AD 
1024 ^3 



256 M 1 ° 



]uf +1) + 0(y 4 ) 



(31) 



That is, to model level £ dynamics at level (•£+!) the coefficients in Q3"0 
transform according to the linear transform 



c («+D =Tc m where T = 



1 
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1 

64 


1 
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^512 


1 
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1 

64 
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16384 


4096 


1034 


256 



(32) 



By induction, the level i dynamics have centred difference coefficients c' £ ' = 
T £ c' ' . Consequently, the dynamics that emerge on macroscale grids are 
determined by the powers T { for large I. Since T is triangular the powers are 
simple: the dominant structure for large I corresponds to that of the leading 
eigenvalue 1 /4; its eigenvector gives the centred difference coefficients that 
emerge on the macroscale as 



AO) 



■l_h __L 1 

M V » 12' 9( 



90' 



1 

560' 



as I 



oo 



(33) 



Recognise in this vector the coefficients of various powers of centred differ- 
ences in a discrete representation of the continuum diffusion operator 3 2 /3x 2 . 
That isj§ continuum diffusion emerges on the macroscale for all isotropic, 
conservative, linear, continuous time, microscale dynamics provided there is 
some component of b 2 in the microscale (c[ ^ ^ 0). 

What is novel here? That continuum diffusion emerges on macroscopic 
scales has been well known for centuries. The novelty is the centre manifold 
theory framework I set up to prove this well known fact. This framework 
illuminates issues and empowers us to analyse much more difficult nonlinear 
dynamics, Section 14.31 and potentially stochastic problems. 

Furthermore, the framework shows that a consistent truncation in the 
interelement coupling parameter y generates a macroscopic approximation 



6 This emergence of diffusion is proved here only up to the order of error in the truncation 
used here. 



21 



to continuum diffusion that is of the same order of error in y. (For exam- 
ple, in the Introduction we discussed multiscale modelling with truncations 
to 0(y 2 ).) This consistency follows because truncating the mapping op- 
erator T, given in ( 1321) . simply truncates its spectrum and truncates its 
eigenvectors (as T is triangular). 

4.2 Renormalise advection-dispersion 

Now consider advection-dispersion on a multiscale hierarchy. Although the 
dynamics are linear, and in contrast to the previous subsection, the transfor- 
mation from one level to another in the hierarchy is nonlinear. 

For example, suppose the microscopic dynamics is simply the discrete 
Uj = [— (0.5 + 5 ]uj . Then computer algebra [32J derives that the multigrid 
hierarchy of dynamic models is 

[-^ + 6 2 ]u|°\ (34) 

\ [ - ^5 + |6 2 + y(-±S 2 + 4^6 3 - ^S 4 )]^ 1 ' + 0(y 2 ) , (35) 

-If — ix5 -h ^6 2 + Yf--6 2 + m u6 3 - 68901 5 4 )W 2) 

4L r 1 ^80 ^rl 80° ^ 3200^° 512000° U U j 

+ 0(y 2 ), (36) 
J[ - m-5 + 0.5001 56 2 

+ y (-0.3751 55 2 + 0.37527^iS 3 - 0.18763S 4 )]uf > + C(y 2 ) ,(37) 
^[-\ib + 0.500008 2 

+ y(-0.498055 2 + 0.49805^6 3 -0.249026 4 )]u| 9) + C(y 2 ) ,(38) 

where for simplicity I omit the level of the discrete mean and difference 
operators. Evidently, as the level I increases, and upon renormalising the 
time scale by the factor of I 1 (the grid step), these models approach a fixed 
point corresponding to an upwind model of the advection. As is well known, 
advection dominates diffusion on large scales. This centre manifold supported 
multiscale transformation preserves the advection speed, and does so with 
stable upwind differencing. 

Now explore the general mapping of linear conservative dynamics from 
one level on the multigrid to another. Generalise the form (1301) for isotropic 



du 



(0) 



dt(°) 

duj 1 ' 

dtd) 

duj 2 ' 

dt« 



du 



(3) 



dt( 3 ) 



du 



(9) 



dt^ 
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dynamics to the general finite difference representation of conservative linear 
operators: 



p=1 I k=1 



(39) 



for some coefficients c Pi k where Cii = — c and C12 = d as used earlier. The 
operator in braces { }, called C v earlier, represents a general operator of stencil 
width 2p + 1 . For example, working to error O (y 3 , y 2 ) , computer algebra [32] 
derives the model at the next level of the the multiscale hierarchy to be, in 
gory detail, 



~dT 
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16c 2 



-c 2)2 + 4c 2 , 4 ) 



Hi 



16c 



1,2 



(— ci.i -2c 2 ,i +4c 2)3 ) + 
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4c 



2,4J 



1,2 



+ o(y 2 ). 



(40) 



This general mapping from (1391) to (1401) governs the particular exapmple hier- 
archy of models (134j) - (138j) . Fine grid scale interactions generate the nonlinear 
dependence upon coefficients shown in the transformation to (HUj) . The ex- 
ample hierarchy (j3*31) - (j3~81) shows that when we scale time by a further factor 
of two in each transformation, to correspond to the time scale of advection of 
a grid of twice the spacing, then the multiscale transformation possess a fixed 
point. Returning to the general transformation from ( 1391) to ( )40l) . rescaling 
time by a factor of two, computer algebra finds precisely two non-trivial fixed 
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points of the multiscale transformation: 
du W 

= c,{^6 + l6 2 + y [-l6 2 ±^-i6 4 ]}uj* ) + (9(y 2 ), 

for some speed c* (positive) which will depend upon the precise microscopic 
system. These fixed points are purely upwind macroscale models of the 
advection and dispersion dynamics. Such stable upwind models naturally 
emerge from our rational transformation of dynamics based upon dynamical 
systems theory. 



4.3 Approximate the nonlinear Burgers' dynamics 



Burgers' partial differential equation, 3u/3t + u 3u/3x = 3 2 u/3x 2 , is fre- 
quently invoked as a benchmark problem in nonlinear spatio-temporal dy- 
namics as it involves the important physical mechanisms of dissipative dif- 
fusion and nonlinear advection. As an example of a nonlinear application of 
our multiscale methodology, suppose Burgers' pde is spatially discretised to 



du 



(0) 



dt(°> 



b 2 



u 



(0) 



(0) c (0) 
ecu UOU- 



(41) 



where the time scale t' ' is chosen to make the coefficient unity for the centred 
difference approximation S 2 Uj of the diffusion 3 2 u/3x 2 . Take equation (141 p 
to be the microscale discrete nonlinear dynamics. The parameter <x measures 
the importance of the nonlinear advection on this microscopic scale. Section [2] 
places the coarse grid modelling of such nonlinear discrete dynamics within 
the purview of centre manifold theory. 

For relatively small parameter a, straightforward modifications of the 
computer algebra for the earlier linear dynamics [32] analyses nonlinear prob- 
lems. The reason is that as long as the nonlinearity is relatively weak, small ct, 
the dominant mechanism in each element is the linear dissipation of S 2 just 
as for the linear dynamics. Our multiscale modelling transforms the fine 
grid dynamics (I4"IT) into the level one dynamics (142]) ; applying the multiscale 
modelling again transforms the level one dynamics (I42p into the level two 
dynamics 



du 



1 1 2 (D 2 

4 + 16 aU ' 
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du 



(2) 



dt« 
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(43) 



Evidently we could continue this transformation across many more levels of 
a multigrid hierarchy 

The transformation from (j4ip to (j42p . to (I43p is based upon small non- 
linearity a. However, as we should expect, the nonlinear advection appears 
to become more important at larger scales: the relative magnitude of the 

■y / f£)\ 2 

nonlinear enhancement to dissipation, a (u- J , increases when going from 
level to level 2. After transforming over enough levels, the nonlinearity 
will begin to dominate the linear basis of the analysis here; at that length 
scale I expect the discrete dynamics to morph into a qualitatively new form, 
one dominated by nonlinear advection. Although these emergent dynam- 
ics cannot be captured by the transformation used here, a generalisation of 
the algebra to being based about a nonlinear subspace of piecewise constant 
solutions may be feasible. In that case, the centre manifold theory of Sec- 
tion [2] would still support the multiscale modelling of the strongly nonlinear 
dynamics. 



5 Conclusion 

This article introduces a new dynamical systems approach to modelling and 
linking dynamics across a multigrid hierarchy. 

Because we recover continuum diffusion, §4.11 and upwind advection, §4.21 
on macroscales we are reassured that the process of modelling from one grid 
to the next coarser grid is indeed sound, as claimed by centre manifold the- 
ory, £J2J Further, errors do not appear to accumulate when we iterate the 
modelling transformation across many changes in length scales. 

At all lengths scales in the hierarchy of models, centre manifold theory 
assures us that the model on each scale is exponentially attractive, §3.31 
and provides an estimate of the rate of attraction. This theoretical support 
applies for the finite spectral gaps on the multigrid hierarchy. 

The geometric picture of invariant slow manifolds also provides a rationale 
for providing initial conditions for the models at each length scale [2U El ETJ, 
e.g.]. Section [33] connects appropriate initial conditions with the restriction 
projection of multigrid solution of linear equations [51 e.g.]. 
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In addition to providing the dynamics at all length scales in the hierarchy, 
this approach also provides intraelement structures realised by the dynamics 
of the grid values: in terms of the level (-6+1 ) dynamic variables, equation ( fT9l 
with ( fiTj) describe the corresponding structures on the level I grid. In some 
sense, equation ( TiTI) provide 'wavelets' for each grid scale P, [161 e -g-]- O ur 
modelling connects the dynamics of wavelets across a hierarchy of length 
scales. 

All the analysis herein is for dynamics in one space dimension. Just as 
for holistic discretisation of pdes [TH], I expect extension to higher space 
dimensions will be straightforward. This article focussed on dynamics which 
to a first approximation could be modelled by advection-dispersion equa- 
tions; extension to dynamics of necessarily higher-order, such as a discrete 
Kuramoto-Sivashinsky equation, could also be analogous to the approach of 
holsitic discretisation [I7[ [18]. Similarly, extension to stochastic mutliscale 
dynamics could be analogous to that of holistic discretisation of stochastic 
dynamics [311 [331 e -g-]- 
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